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The van der Waals quintessence equation of state is an interesting scenario for describing the 
late universe, and seems to provide a solution to the puzzle of dark energy, without the presence 
of exotic fluids or modifications of the Friedmann equations. In this work, the construction of 
inhomogeneous compact spheres supported by a van der Waals equation of state is explored. These 
relativistic stellar configurations shall be denoted as van der Waals quintessence stars. Despite 
of the fact that, in a cosmological context, the van der Waals fluid is considered homogeneous, 
inhomogeneities may arise through gravitational instabilities. Thus, these solutions may possibly 
originate from density fluctuations in the cosmological background. Two specific classes of solutions, 
namely, gravastars and traversable wormholes are analyzed. Exact solutions are found, and their 
f— ^ respective characteristics and physical properties are further explored. 
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| The universe is presently undergoing an accelerated phase of expansion. Several candidates have been proposed in 
. the literature to explain this cosmic accelerated expansion, such as dark energy models, the generalized Chaplygin 
gas, modified gravity and scalar-tensor theories, tachyon scalar fields and specific braneworld models, such as the 
Dvali-Gabadadze-Porrati (DGP) model. However, it was pointed out that assuming that dark energy is governed by 
a perfect fluid equation of state may systematically induce wrong results and be misleading in inferring the nature of 
dark energy In this spirit, an alternative model was recently proposed without the presence of exotic fluids and 
modifications of the Friedmann equations, using a more complicated equation of state, namely, the van der Waals 
( VDW) equation of state 0, Q , given by 
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& p = T^T p - ap > (1) 

bJ[). where p is the energy density and p the pressure of the VDW fluid. The accelerated and decelerated periods depend 
on the parameters, a, f3 and 7 of the equation of state, and in the limiting case a, (3 — ► 0, one recovers the dark energy 
• i-H . equation of state, with 7 = p/p < — 1/3. It was also stressed that the perfect fluid equation of state p — reflects 
an approximation of cosmic epochs describing stationary situations, where phase transitions are not considered (4j. 
Thus, one of the advantages of the VDW model is that it describes the transition from a matter field dominated era 
to a scalar field dominated epoch, without introducing scalar fields. Furthermore, it is useful to explain the universe 
with a minimal number of ingredients, and the VDW gas actually treats dark matter and dark energy as a single fluid. 
The VDW quintessence scenario has also been confronted successfully with a wide variety of observational tests, by 
constraining the free parameters [1] . 

The success of this model has stimulated several different approaches. Recently, a model with a binary mixture 
whose constituents are described by a VDW fluid and by a dark energy density was also proposed [5(] . It is interesting 
to note that this model can simulate several aspects, namely, an inflationary period where the acceleration grows 
exponentially; and a present accelerated period where the dark energy density dominates over the energy density of 
the VDW fluid. The construction of a general scheme where the VDW ideal fluid has a mathematically equivalent 
representation as a scalar-tensor theory with a specific potential, as well as explicit examples were also explored in 
Ref. (f|. Another approach has been analyzed in the context of brane cosmology pj, where the cosmological fluid 
on the brane is modeled by the VDW equation of state. It was shown that this model reproduces several features, 
namely, an initial accelerated epoch where the VDW fluid behaves like a scalar field with a negative pressure, and a 
present accelerated phase due to a cosmological constant on the brane. 
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Despite of the fact that the VDW equation of state represents a spatially homogeneous cosmic fluid and is assumed 
not to cluster, due to gravitational instabilities inhomogeneities may arise. Thus, the solutions outlined in this paper 
may possibly originate from density fluctuations in the cosmological background, resulting in the nucleation through 
the respective density perturbations. Therefore, the pressure in the VDW equation of state may be regarded a radial 
pressure, and the transverse pressure is determined through the Einstein field equations. The justification of this 
approach may be considered as an extension of the construction of inhomogeneous wormhole spacetimes supported 
by phantom energy 0. The latter construction was motivated by the analysis carried out in Ref. 0], where time- 
dependent solutions were constructed with ghost scalar fields. In this model, it was shown that the radial pressure is 
negative everywhere and equals the transverse pressure far from the throat, showing that ghost scalar fields behave 
essentially as dark energy. Now, in the VDW context, one may also consider, in principle, that for large radial distances 
in these inhomogeneous solutions the radial pressure equals the transverse pressure, and thus behaving essentially as a 
VDW cosmological fluid, consequently justifying the above approach. The transverse pressure obtained from the field 
equations should not be considered typical of the cosmic fluid driving the accelerated expansion, as it is only a feature 
of these inhomogeneous VDW configurations. It is also important to emphasize that the equation of state leading to 
the acceleration of the Universe on large scales, is an average equation of state corresponding to a background fluid. So 
in the inhomogeneities giving rise to the configurations discussed in this work, it is also important to justify if one may 
really use an equation of state that in fact corresponds to a background fluid. For this, it is important to understand 
the evolution of structure formation in the VDW scheme. It was shown that for large redshifts the VDW model is 
similar to the ACDM model [J], where the VDW dark matter plays the role of a cosmological constant, showing that 
structure formation could evolve in a similar manner as in the ACDM model, although these points need to be further 
investigated Q. Now, it was also shown that due to a negative adiabatic sound speed, perturbations are unstable 
in the VDW quintessence scenario, and inevitably the role of anisotropic stresses and entropy production has to be 
taken into account. Thus, as a first approximation and as a working hypothesis, one may argue that this also justifies 
the above approach of using an equation of state that corresponds to a background fluid, with a radial pressure. 
We shall denote these inhomogeneous spherically symmetric solutions as van der Waals quintessence stars. It is also 
interesting to note that the thermodynamic properties of a general relativistic isotropic pressure self-gravitating ball 
of fluid, undergoing a generic first order phase transition and governed by the VDW equation of state, was studied in 
Ref. [10( • In this work, we shall generalize the latter to an anisotropic pressure general relativistic model, and we shall 
consider two specific class of solutions of VDW quintessence stars, namely, gravastars and traversable wormholes. 

The "gravastar" (gravitational vacuum star) picture developed by Mazur and Mottola [llj , initially arouse from a 
considerable amount of scepticism related to spacetime singularities and event horizons in the Schwarzschild black 
hole solution (see Ref. [l2j for a review on specific misconceptions and ambiguities related to these issues). In the 
gravastar model, the interior Schwarzschild solution is replaced with a de Sitter condensate, and thus, does away with 
the problem of the singularity at the origin and the event horizon. In this model, the quantum vacuum undergoes a 
phase transition at or near the location where the event horizon is expected to form. In this spirit, it is interesting to 
note that several models have been proposed to some extent in the literature [l3[ . Recently, motivated by dark energy 
quintessence, which is a possible candidate responsible for the late-time cosmic accelerated expansion, a generalization 
of the gravastar model has been proposed [l5| , by considering an interior solution governed by the dark energy 
equation of state, 7 = p/p < —1/3. Stable dark energy stellar models were analyzed in Ref. 15], and it was found 
that large stability regions exist that are sufficiently close to where the event horizon is expected to form. It was 
also emphasized that it would be difficult to distinguish the exterior geometry of the dark energy stars from an 
astrophysical black hole. In this context, it is interesting to note that stars supported by a generalized Chaplygin gas 
equation of state have been explored [la ]. 

In this work, we shall also be interested in investigating the construction of traversable wormholes [171 ] supported 
by the VDW equation of state. An analogous approach was recently carried out in the presence of phantom energy 
3 and the generalized Chaplygin gas [la] . We will show that traversable wormhole solutions may be constructed 
using the VDW equation of state, which are either asymptotically flat or possess finite dimensions, where the exotic 
matter is confined to the throat neighborhood. The latter solutions are constructed by matching an interior wormhole 
geometry to an exterior Schwarzschild vacuum. Analogously to their phantom and Chaplygin counterparts, these 
VDW quintessence wormholes have far-reaching astrophysical and cosmological consequences, such as the production 
of closed timelike curves and the consequent violation of causality. 

This paper is outlined in the following manner: In Sec. the equations of structure for van der Waals quintessence 
stellar models are presented. In Sec. [31 an exact solution of a gravastar is found, and the respective properties and 
characteristics are analysed. In Sec. [H an exact solution of a traversable wormhole geometry is also found. In Sec. (5[ 
we conclude. 



3 



2. VDW QUINTESSENCE STARS: EQUATIONS OF STRUCTURE 

Consider a static and spherically symmetric spacetime, given by the following metric, in curvature coordinates 

ds 2 = -e 2 *^ dt 2 + ^-y-- + r 2 (d9 2 + sin 2 9 dej?) , (2) 

1 — 2m\r)/r 

where <£>(r) and m(r) are arbitrary functions of the radial coordinate, r. The function m(r) is the quasi-local mass, 
and is denoted as the mass function. 

The Einstein field equation, G^ v = 87rT^„, where is the Einstein tensor and T^ v the stress-energy tensor, 
provides the following relationships 

m = 4:Trr 2 p , (3) 

*'=^±^, (4) 
r(r — 2m) 

(p+p r )(m + 47rr 3 p r ) 2 f 

Pr = 7 7r~\ + -(Pt-Pr), 5 

r(r — 2m) r 

where the prime denotes a derivative with respect to the radial coordinate, r. p(r) is the energy density, p r (r) is the 
radial pressure, and pt(r) is the lateral pressure measured in the orthogonal direction to the radial direction. Equation 
corresponds to the anisotropic pressure Tolman-Oppenheimer-Volkoff (TOV) equation, and may be obtained using 
the conservation of the stress-energy tensor, T^ v . v = 0. 

The four-velocity of a static observer, at rest at constant r,9,<p, is = dx^/dr = ([/', 0,0,0) = (e~*( r ), 0, 0, 0). 
The observer's four-acceleration is a M = [/ M ; „ U v , so that taking into account metric ((2]), we have a* = and 

a r = T\ t (I-) 2 =*'(l-6/r). (6) 

Note that from the geodesic equation, a radially moving test particle, which initially starts at rest, obeys the following 
equation of motion 

d 2 r „, f dt 



d r*-- Tr "{Tr) =~ ar - W 

a r is the radial component of proper acceleration that an observer must maintain in order to remain at rest at constant 
r, 8, cf). One may consider that the geometry is attractive if a r > 0, i.e., observers must maintain an outward-directed 
radial acceleration to keep from being pulled into the star; and repulsive if a r < 0, i.e., observers must maintain an 
inward-directed radial acceleration to avoid being pushed away from the star. This distinction depends on the sign of 
<&', as is transparent from Eq. ©. In particular, for a constant redshift function, <J>'(r) = 0, static observers are also 
geodesic. Thus, the convention used is that $'(r) is positive for an inwardly gravitational attraction, and negative for 
an outward gravitational repulsion (l9j |. 

Despite of the fact that, in a cosmological context, the van der Waals fluid is considered homogeneous, inho- 
mogeneities may arise through gravitational instabilities, resulting in a nucleation of the cosmic fluid due to the 
respective density perturbations. Thus, these solutions may possibly originate from density fluctuations in the cos- 
mological background. Consider the VDW equation of state for an inhomogeneous spherically symmetric spacetime, 
given by 

IP 2 , al 

Pr = T^Tp ap ■ (8) 

In the present context of relativistic stellar models, we shall consider the pressure in Eq. (JSJ) as the radial pressure, 
p r , and the tangential pressure may be obtained via Eq. (J5]). Note that the limiting case a, (3 — > 0, reduces to the 
dark energy equation of state, with 7 < —1/3. Using Eqs. ©-(HI), then Eq. ([8]) provides the following relationship 

2m\ , m 7m' am' 2 , , 

r r 1 - 2™' 4 7rr 2 v 1 

Note that now we have four equations, Eqs. ©-(O and Eq. ([5]), with five unknown functions of r, i.e., 3>(r), m(r), 
p(r), p r {r) and pt{r). To solve the system, one may adopt different strategies, namely, one may model an appropriate 
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spacetime geometry by imposing m(r) and/or $(r) by hand and consequently determine the stress-energy tensor 
components. In counterpart, one may consider an adequate source of the spacetime geometry by imposing the stress- 
energy components, and consequently determine the metric fields. 

Another useful and interesting relationship may be obtained from Eqs. ([3]), ([5]) and Eq. (jHJ), which provide 
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(10) 



A = p t — p r is denoted the anisotropy factor, as it is a measure of the pressure anisotropy of the fluid comprising the 
VDW quintessence star. The factor A/r represents a force due to the anisotropic nature of the stellar model, which 
is repulsive, i.e., being outward directed if pt > p r , and attractive if pt < p r . A = corresponds to the particular case 
of an isotropic pressure VDW quintessence star. 

One may, in principle, construct asymptotically flat spacetimes, where $(r) — * 0, and m(r)/r — > 0, as r — > oo. An 
alternative approach is to consider a cut-off of the stress-energy tensor at a junction radius a. For instance, consider 
for simplicity that the exterior solution is the Schwarzschild spacetime, given by 



ds 2 



2M\ , 9 ( 2M , , ., ., 
1 - — ) clt 2 + ( 1 - — ) dr 2 + r 2 (d0 2 +sin 2 6» ( 



(11) 



M may be interpreted as the VDW quintessence star's total mass. In this case the spacetimes given by the metrics 
Eq. and (jlip are matched at a, and one has a thin shell surrounding the star. Using the Darmois- Israel formalism, 
the surface stresses are given by 

/ / 9 ]\,f / r,™l„\ \ 

(12) 

J) _ ^ a ' " _ J — " - \u-^ni) (13) 

8na 

where the overdot denotes a derivative with respect to the proper time, r; a is the surface energy density and V the 
surface pressure (see Refs. [2(| for details). The static case is given by taking into account a = d = 0. The total mass 
of the VDW quintessence star, for the static case, is given by 




M = m(a ) + m s (a ) 



2m(a ) _ m s (a ) 
a 2ao 



(14) 



where 



is the surface mass of the thin shell, and is defined as to, = Aira 2 a. 



3. VDW QUINTESSENCE GRAVASTARS 



The Mazur-Mottola model is constituted by an onion-like structure with five layers, including two thin-shells, with 
surface stresses a± and V±, where a is the surface energy density and V is the surface tangential pressure. The 
interior of the solution is replaced by a segment of de Sitter space, which is then matched to a finite thickness shell 
of stiff matter with the equation of state p = p. The latter is further matched to an external Schwarzschild vacuum 
with p = p = 0. A simplification of the Mazur-Mottola configuration, and its dynamic stability, was considered in 
Ref. [21( using a three-layer solution [2l| . i.e., a de Sitter interior solution was matched to a Schwarzschild exterior 
solution at a junction surface, comprising of a thin shell with surface stresses a and V '. 

Consider the specific case of a constant energy density, p(r) = po, so that Eq. ([3]) provides the following mass 
function 

m(r) = Ar 3 , (15) 
where for simplicity, the definition A = 4irpo/3 is used. Using Eq. ©, one finds the following expression 
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where, for notational convenience, the constant \ is defined by 
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Now, in order to obtain an accelerated behaviour of the Universe, the condition p + 3p < 0, should be obeyed. We 
shall use this condition, in the present inhomogeneous spacetime, which yields the following inequality 
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(18) 



Note that in order to provide a positive energy density, one may also impose the following simultaneous conditions, 
in terms of the VDW equation of state parameters 
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See Ref. Q for further details. Taking into account the mass function given by Eq. (jT5j) , and using Eq. ([3|), we have 
p = p = 3 A/ (Air), so that inequality (fT5)) takes the following form 

3A X 
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< 0. 



(20) 



from which one readily verifies, using Eq. (|16|) . that $'(r) < 0, providing the necessary repulsive character, which is a 
fundamental property of gravastar models. Note that as we are considering positive energy densities, so that A > 0, 
in the gravastar models, we have \ < 0. 

Equation (|T6)) . can be integrated to provide the following spacetime metric 



ds 2 = - (1 - 2A 
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The stress-energy scenario is given by p = po = 3A/(A-k) and 
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The anisotropy factor is provided by the following relationship 

A 2 r 2 X {2 + X ) 
8tt(1 - 2Ar 2 ) ' 

which is plotted in Fig. [TJ Note that A = 0, i.e., p r — p t , at the center r — as was to be expected. 
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4. VDW QUINTESSENCE TRAVERSABLE WORMHOLES 



The spacetime metric representing a spherically symmetric and static wormhole is given by Eq. l[2]). where <&(r) 
and b(r) = 2m(r) are arbitrary functions of the radial coordinate, r, denoted as the redshift function, and the form 
function, respectively (Tt| . The radial coordinate has a range that increases from a minimum value at ro, corresponding 
to the wormhole throat, to infinity. To be a wormhole solution a flaring out condition of the throat is imposed, i.e., 
(b — b'r)/b 2 > [rr| . From this we verify that at the throat b(r ) = r = r , the condition b'(r ) < 1 is imposed to 
have wormhole solutions. For the wormhole to be traversable, one must demand that there are no horizons present, 
which are identified as the surfaces with e 2 * — > 0, so that $(r) must be finite everywhere. Note that the condition 
1 — 6/V>0is also imposed. 

Consider, for instance, the specific case of of b(r) — rfi/r, which obeys the conditions of a wormhole, so that Eq. 
© provides the following solution 
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FIG. 1: Plot of the anisotropy factor for a VDW quintessence star with a constant energy density. We have defined the following 
dimensionless parameters: S — A/ A and a = \J~Ar. 



where C is a constant of integration. We verify the existence of an event horizon, r = ro, 
right-hand-side. Thus, fine-tuning the 7 parameter as 
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This now corresponds to a traversable wormhole, as one verifies the absence of an event horizon. Note that this 
solution is not asymptotically flat, so that one needs to match the latter to an exterior vacuum spacetime. If we 



consider a further imposition, namely, a = am; 
assumes a particularly simple form, given by 



0' 



then Eq. 



becomes zero, 7 = 0, so that the redshift function 
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The VDW equation of state, Eq. {2} , for this particular choice of parameters is p r - 
radial pressure. Note that this specific choice of a wormhole is asymptotically flat, 
take the following form 



(28) 

— 87rrQ p 2 , implying a negative 
The stress-energy components 
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(29) 
(30) 



This solution corresponds to a rather exotic form of a VDW fluid, due to the presence of a negative energy density. 
The null energy condition, defined as T^ v k u > 0, where fc M any null vector, is also violated, i.e, p + p r < 0, which 
is a necessary condition of wormholes. 

It is also of interest to note that the VDW cosmic fluid with an isotropic pressure may be completely determined 
by only one parameter in the equation of state, so that one may consider 7 as the only independent parameter needed 
to describe the VDW fluid [|[. Therefore, using a wide variety of values of 7 constrained by observations considered 
in Ref. @, one may obtain restrictions on the parameters a and j3, through Eq. (|26p . in the inhomogeneous solutions 
considered here. On the other hand, although the particular case of 7 = is not included in the best fits parameter 
range Q, it is not excluded from the observational constraints. This should not be considered a serious shortcoming, 
as the parameters a, (3 and 7 are defined using the critical values of the energy density, the isotropic pressure and 
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the critical volume [|| , and in the presence of anisotropic pressures the respective definitions should be extended to 
include the transverse pressure. Despite of this fact, and due to the analytical complexity of Eq. (f2"T)) . we shall, 
for simplicity, consider the case 7 = and further analyze the physical properties and characteristics of this specific 
wormhole solution below, which prove to be extremely interesting. 

One may also quantify the "total amount" of energy condition violating matter, which amounts to calculating the 
"volume integral quantifier" defined as J T^ u k^k v dV [22j|. The amount of violation is defined as the extent to which 
this integral becomes negative. Using the "volume integral quantifier" , which provides information about the "total 
amount" of averaged null energy condition (ANEC) violating matter in the spacetime (see Ref. [22j for details), given 
by Iy = J [p(r) +p r (r)] dV, with a cut-off of the stress-energy at a, we have 



Iv = / (r-6) 



In 



1 -b/r 



dr. (31) 



Now, using the form function given by b(r) — r^/r, and the redshift function provided by Eq. (|28D . and evaluating 
the integral, one finally ends up with the following simplified expression for the "volume integral quantifier" 

r (r% + 5r a 4 - 6a 5 ) 
Iy - 5a^ ■ (32) 

By taking the limit a — > r , one readily verifies that Iy — > 0. It is also interesting to note that in the limit a — > oo, 
the volume integral tends to a finite value, i.e., Iy = — 6ro/5. This proves that as in the specific case of phantom 
[| and Chaply gin [l8j wormholes, one may theoretically construct a wormhole with arbitrarily small amounts of a 
VDW quintessence fluid. As emphasized in Ref. [8|, this result is not unexpected, however, it is interesting to note 
the relative ease with which one may theoretically construct wormholes supported by infinitesimal amounts of exotic 
fluids used in cosmology to explain the present accelerated cosmic expansion. 

It is also of a particular interest to analyze the traversability conditions. We will be interested in specific solutions 
for traversable wormholes and assume that a traveller of an absurdly advanced civilization, with human traits, begins 
the trip in a space station in the lower universe, at proper distance I — —h, and ends up in the upper universe, at 
I = The proper distance is given by dl — ±(1 — bjr)^ 1 ! 2 dr. Assume that the traveller has a radial velocity v(r), 
as measured by a static observer positioned at r. One may relate the proper distance travelled dl, radius travelled dr, 
coordinate time lapse dt, and proper time lapse as measured by the observer dr, by the following relationships 

dl ( _ b\~ 1/2 dr 



s = Te [ 1 -- r > *• < 33 > 

, 7 * ^-iV 1 "* ,34, 
dr \ r J dr 

See [TtJ for details. 

For a convenient trip through the wormhole, certain conditions should also be imposed [T^. Firstly, the entire 
journey should be done in a relatively short time as measured both by the traveller, Ar tr , and by observers who 
remain at rest at the stations, At st . For simplicity, we shall take into account non-relativistic and constant traversal 
velocities, 7 w 1. Thus, Art r and At s t are given by 

A _ r _ f +h dl 2 r dr >2 f° dr w 2(a - rg) ^ 

J-h «7 Jro (1 - b/r) 1 / 2 vj ~ J ro vj v 
f +h dl A" dr [ a dr 2 f a a 

Considering the redshift function provided by Eq. (|28|) . and a — 2r$, so that dr w 1.18 ro, we verify that both 

Artr and At B t are bounded from below by ~ 2r§jv. 

An important traversability condition required is that the acceleration felt by the traveller should not exceed Earth's 
gravity [171 ] . Thus, the traveller's four-acceleration expressed in his proper reference frame, considering once again 
non-relativistic and constant traversal velocities, 7 ~ 1, yields the following restriction (see Ref. [ItJ for details) 



b\ 1/2 , 
1 - - 1 

r 



< ffffi , (37) 



which is readily verified at the throat, where the conditions are most severe. 
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Another important condition is that an observer traversing through the wormhole should not be ripped apart 
by enormous tidal forces. Thus, it is required that the tidal accelerations felt by the traveller should not exceed, 
for instance, the Earth's gravitational acceleration [r?| . The constraint |Aa^ | < provides the tidal acceleration 
restrictions as measured by a traveller moving radially through the wormhole, given by the following inequalities 



l- b - 

r 
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2r 2 



($')2 _ 



b'- - 



b'r-b 
2r(r - b) 



2(r-6)$' 



< 



3© i 



W I < 5© • 



(38) 
(39) 



We refer the reader to Rcf. 

.17] for details related to the deduction of these conditions. The factors j^ 1 | and \rf \ 
appearing in inequalities (f3"5| - (|3"9"|) are the spatial separations between radial and lateral parts of the traveller's body, 
as measured in his proper reference frame. For computational purposes we may assume 1 77 1 | w |?7 2 | ps 2 m. 

The radial tidal constraint, Eq. (f3"5| . constrains the redshift function, and the lateral tidal constraint, Eq. (|39[) . 
constrains the velocity with which observers traverse the wormhole. These inequalities are particularly simple at the 
throat, ro, and reduce to 



l*'(ro)| < 



7 D < 



2ff© r 



a- wr 

2ff© rj 

{i-v)W'\- 



(40) 
(41) 



Taking into account the redshift function given by Eq. (|2"8"]l . and the form function considered above, i.e., b(r) = r^/r, 
we verify that evaluated at the throat we have <&'(ro) = 1/ro and fr'^o) = —1. Thus, inequality (|40[) provides the 
restriction 1 < ge^/l 7 ? 1 I- Considering the minimum value, i.e., <?e r o/l 7 ? 1 I ~ 1; an d that {rj 1 | ss \q 2 | sa 2m, we 
verify that the minimum value of the throat is given by ro ~ 2 x 10 8 m. Using the latter approximation, and taking 
into account non-relativistic and constant traversal velocities, 7 « 1, then inequality (|4ip provides the following 
restriction for the velocity, v < c. For instance, considering that the traversal velocity be of the order v ~ 10 -2 c, 



and assuming that the traversal times take their minimum value, we have Ar tr ~ 
extremely reasonable traversable wormhole, useful for travellers with human traits. 



10 s. This provides an 



5. CONCLUSION 



In conclusion, noting that the VDW quintessence equation of state is an interesting scenario for describing the 
late universe, we have explored the construction of relativistic stellar models, in particular gravastars and traversable 
wormholes, supported by a VDW equation of state. These solutions may possibly originate from density fluctuations 
in the cosmological background, resulting in the nucleation through the respective density perturbations. We have 
found an exact solution of a gravastar model, by considering the specific case of a constant energy density. Relatively 
to traversable wormhole geometries, we also found an exact solution of an asymptotically flat wormhole, as well as 
solutions, where the exotic matter is constrained to the throat neighborhood, by considering a matching of an interior 
wormhole geometry to an exterior Schwarzschild vacuum solution. It was also found that one may theoretically 
construct a wormhole with arbitrarily small amounts of a VDW quintessence fluid. Several characteristics and 
properties of these traversable wormholes using the traversability conditions were also explored. 
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